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Abstract 

We study the dynamical chaos and integrable motion in the planar circular restricted three-body problem and determine the fractal 
dimension of the spiral strange repeller set of non-escaping orbits at different values of mass ratio of binary bodies and of Jacobi 
'O integral of motion. We find that the spiral fractal structure of the Poincare section leads to a spiral density distribution of particles 
remaining in the system. We also show that the initial exponential drop of survival probability with time is followed by the algebraic 
decay related to the universal algebraic statistics of Poincare recurrences in generic symplectic maps. 

^ Keywords: Chaos, Celestial mechanics, Binaries, Galaxies: spiral 


1. Introduction 

I The restricted three-body problem was at the center of stud- 
CIh ies of dynamics in astronomy starting from the works of|Euler| 


( |1772| ), pacobi| ( |1836| ) and |Poincare| ( |1890| ). The progress in the 
cH understanding of this complex problem in XXth and XXIth cen- 
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Q^turies is described in the fundamental books { Szebehely 1967t 

^ |Henon[ |1997[ |2001t [Valtonen and Karttunen[ 2006| ). As it was 

proven by Poincare] ( |1890| ) in the general case this system is 
not integrable and only the Jacobi integral is preserved by the 
dynamics ( jJacobij |1836| ). Thus a general type of orbits has a 
chaotic dynamics with a divided phase space where islands of 
stability are embedded in a chaotic sea (Chirikov 1979[ Licht- 


(N 

> 

00 

m 

i> 

o 

o\ 

o 


jenberg and Lieberman||1992t|Qtt||1993|). 

In this work we consider the Planar Circular Restricted 
Three-Body Problem (PCRTBP). This is an example of a con¬ 
servative Hamiltonian system (in a synodic or rotating reference 
frame of two binaries) with two degrees of freedom. However, 
this is an open system since some trajectories can escape to 
infinity (be ionized) so that the general theories of leaking sys¬ 


tems (Altmann et al. 20131 and naturally open systems (e.g. 


jContopoulos and Efstathiou 


|2004| ) find here their direct appli- 




cations. It is known that such open systems are characterized 
by strange repellers related to non-escaping orbits and by an 
exponential time decay of probability to stay inside the system. 
However, as we show, in the PCRTBP system with a divided 
phase space one generally finds an algebraic decay of proba¬ 
bility of stay related to an algebraic statistics of Poincare re¬ 
currences in Hamiltonian systems (see e.g. jChirikov and She- 
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effect appears due to long sticking of trajectories in a vicin¬ 
ity of stability islands and critical Kolmogorov-Amold-Moser 
(KAM) curves. Thus an interplay of fractal structures and alge¬ 
braic decay in the PCRTBP deserves detailed studies. 

Among the recent studies of the PCRTBP we point out the 
advanced results of jNagler ( |2004| |2005| ) where the crash prob¬ 
ability dependence on the size of large bodies has been studied 
and the fractal structure of non-escaping orbits has been seen 
even if the fractal dimensions were not determined. This re¬ 
search line was extended in Astakhov and Earrelly ( 2004| ); As- 
takhov et~ST] ( |2005| ) with a discussion of possible applications 
to the Kuiper-belt and analysis of various types of orbits in jBar-j 
rio et al.| ( |2009| ); [^tos| ( |2015| ). The analysis of orbits in three 
dimensional case is reported in jMako et al.| ( |2Q10| ) and basin of 
escaping orbits around the Moon has been determined in de As- 
|sis and Te^^ ( |2014| ). 

In this work we determine the fractal dimension of non¬ 
escaping orbits for the PCRTBP with comparable masses of 
heavy bodies and consider the properties of Poincare recur¬ 
rences and the decay probability of stay in this system. The 
system description is given in Section 2, the structure of strange 
repeller is analyzed in Section 3, the decay of Poincare recur¬ 
rences and probability of stay are studied in Section 4, a sym¬ 
plectic map description of the dynamics is given in Section 5, 
discussion of the results is presented in Section 6. 


2. System description 

The PCRTBP system is composed of a test particle evolving 
in the plane of a circular binary whose primaries have masses 
mi = 1 - yu and m 2 = with mi > m 2 . In the synodic frame the 
dynamics of the test particle is given by the Hamiltonian 


H(px, Py 


. ^ {pI + pf) + ypx - xpy + V(x,y) (1) 
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Figure 1: (x, y) - Poincare sections of the Jacobi Hamiltonian ^ with f = 0 and 
0 < 0. Poincare sections for a binary with jj, = 0.3, C = 3 are shown in panels: 
(a) at a large scale, ib) at an intermediate scale, (c) close-up in the vicinity of 
the primary mass. Panel (d) shows the Poincare section for ju = 0.5 and C = 3. 
Red regions are forbidden since there x^ + < 0. Black dots represent non- 

escaped orbits staying inside the r < = 10 region after time r = 10. Invariant 

KAM curves (blue dots) are obtained choosing initial conditions inside KAM 
islands. The red (blue) star * (*) gives the position of the 1 - ju mass (// mass). 
The Poincare section is obtained with orbits of iV = 10^ test particles initially 
placed at random in the region 1.3 < r < 2.5. Particles as escaped once r > Rs. 


where x and y are the test particle coordinates, = x - y and 
Py = y + A are the corresponding canonically conjugated mo¬ 
menta, and 


V(x,y) = - 


(l-yu) 




1/2 


( 2 ) 


is the gravitational potential of the two primaries. Here the dis¬ 
tance between primaries is 1, the total mass mi + m 2 = 1, the 
gravitational constant Q = consequently the rotation period 
of the binary is 2n. Hamiltonian ^ with p otential ^ repre¬ 
sents the Jacobi integral of motion ^acobil 1836 ). In the fol¬ 
lowing we define the Jacobi constant as C = -2H. This Jacobi 
Hamiltonian describes also the planar dynamics of an electri¬ 
cally charged test particle experiencing a perpendicular mag¬ 
netic field and a classical hydrogen-like atom with a Coulomb¬ 
like potential 

We aim to study the dynamics of particles evolving on es¬ 
caping and non-escaping orbits around the binary. We perform 
intensive numerical integration of the equations of motion de¬ 
rived from Hamiltonian Q using an adaptive time step 4th or¬ 
der Runge-Kutta algorithm with Levi-Civita regularization in 
the vicinity of the primaries ( |Levi-Ci\^ 1920| ). The achieved 
accuracy is such as the integral of motion relative error is less 
than 10“^ (10“^) for more than 91% (99%) of integration steps. 
For different Jacobi constants C, we randomly inject up to 10^ 


X X 

Figure 2: (px, x) - Poincare section of the Jacobi Hamiltonian withy = 0 and 
Py > 0 for a binary with p = 0.3, C = 3 (corresponding to Fig. \a, b, c). Panel 
(a): Poincare section at large scale; panel (b): zoom in the vicinity of primaries. 
Black dots represent non-escaped orbits staying inside the r < = 10 region 

after time t = 10. Blue dots represent bounded orbits inside stability islands. 
The red (blue) star * (*) gives the position of the primary (secondary) mass as 
in Fig. 1. The Poincare section is obtained with the same orbits as in Fig. 1. 


test particles in the 1.3 < r < 2.5 ring with initial radial veloc¬ 
ity r = 0 and initial angular velocity 0 < 0 (r and (p are polar 
coordinates in the synodic frame). Each test particle trajectory 
is followed until the integration time attains tg = 10"^ or until 
the region r > Rg = 10 is reached where we consider that test 
particles are escaped (ionized) from the binary. 


3. Strange repeller structures 


In phase space, orbits are embedded in a three-dimensional 
surface defined by the Jacobi constant C. In order to moni¬ 
tor particle trajectories we choose a two-dimensional surface 
defined by an additional condition. Here we choose either the 
condition (r = 0,0 < 0) to represent Poincare section as a (v, y)- 
plane (Figs. 1,5,6,7 and 10) or the condition (y = Q^py > 0) to 
represent Poincare section as a {p^, v)-plane (Fig. 2). A similar 


approach was also used in |Nagler| ( |2004||2005] ). 

We show in Fig. 1 (panels a,b,c) an example of (v,y) - 
Poincare section of the Jacobi Hamiltonian Q for mass param¬ 
eter p = 03 and Jacobi constant C = 3. Red regions correspond 
to forbidden zones where particles would have imaginary ve¬ 
locities. Inside central islands in the close vicinity of primaries 
blue points mark out regular and chaotic orbits of bounded mo¬ 
tion. In particular, the KAM invariant curves ( [Lichtenberg and 
Lieberman| |1992| ) can be seen e.g. in Fig. Ic. In Fig. la, the 
trace of non-escaped chaotic orbits (black points) remaining in¬ 
side the disk r < Rg = 10 after time t = 10 defines a set of 
points forming two spiral arms centered on the binary center of 
mass. This set has a spiral structure of strange repeller since 
orbits in its close vicinity rapidly move away from the set and 
consequently from the binary. The fractal property results in 
a self-similar structure clearly seen by zooming to smaller and 
smaller scales (see Fig. la, b, c). This fractal structure remains 
stable in respect to moderate variation of mass ratio p as it is 
illustrated in Fig. 1 (panels b,d). A strange repeller structure 
is also clearly visible in the corresponding Poincare section in 
{px, x) plane shown in Fig. 2 for yu = 0.3, C = 3. 
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Figure 3: Number of boxes Ni, covering at scale b non-escaped orbits structure 
(strange repeller) appearing in (x, y) - Poincare section of the Jacobi Hamilto¬ 
nian Box-counting computation is performed in an annulus square consist¬ 
ing in a square region -3.9 < x,}; < 3.9 deprived of its central square region 
-1.3 < x,y < 1.3. The annulus square is divided into 8 equal square areas 
of linear size dlo = 1.3. We average the box counting Nt over the 8 squares. 
At scale b each square is divided into Ijb^ boxes of linear size dl = bdlo. 
Box-counting results are shown for Poincare section of orbits staying in the 
r < Rs = 10 disk after t = 3 (black crosses), t = 10 (red squares), = 30 (green 
circles), t = 50 (blue triangles). The fractal dimension D of the strange repeller 
is the slope of the affine function log 2 = /(logiCV^))- Keeping orbits stay¬ 
ing in the r < Rs disk after time ? = 10 (/^ = 3) we obtain a strange repeller 
fractal dimension (a) D = 1.8711 ± 0.0100 (D = 1.8732 + 0.0105) for fi = 0.3 
and C = 3 (see Fig. la,b,c), (b) D = 1.8657 + 0.0117 (D = 1.8690 + 0.0129) 
for // = 0.4 and C = 3 (see Fig. lOb), (c) D = 1.8700 + 0.0077 (D = 
1.8722 + 0.0084) for // = 0.3 and C = 2.6 (see Fig. 5a), (d) D = 1.8349 + 0.0484 
(D = 1.8464 + 0.0436) for ju = 0.3 and C = 3.4 (see Fig. 5d). Fits have been 
performed in the scale range 2^ < Ijb < 2^. We used A = 10^ (a,b,c), 
N = 10^ (d) test particles initially distributed at random in the 1.3 < r < 2.5 
ring. The fractal dimension has been computed with (a) Nt >3 = 39526570, 
Nt>io = 9933333, Nt >30 = 768282, Nt> 5 o = 83290 points, (b) Nt >3 = 
26743797, A,>io = 6550163, A, >30 = 372871, A, >50 = 25037 points, (c) 
Nt >3 = 37610948, A^>io = 8721338, Nt >30 = 419296, Nt> 5 o = 39891 points, 
(d) Nt >3 = 8569720, A, >10 = 5447406, A, >30 = 2245927, A, >50 = 1083887 
points. The curves bend down for the smallest scales b when Ijb^ becomes of 
the order of the number of collected points in Poincare section. 


The fractal dimension of the strange repeller is determined 
using box-counting method (|Lichtenberg and Lieberman||1992[ 
|Ottl[T99^ as Z) = lim/^^o In In (1 /b) where is at scale b 
the number of at least once visited boxes in the Poincare sec¬ 
tion. The box-counting fractal dimension of the strange repeller 
presented in Fig. \a,b,c D ^ 1.87 (Fig. ?>a). This fractal 
dimension value is computed from the strange repeller struc¬ 
ture formed by orbits staying in the disk r < Rs = 10 after 
time t = 10. We see from Fig. 3 that the fractal dimension re¬ 
mains practically the same for the parameters considered here 
A = 0.3,0.4 and C = 2.6,3,3.4. When the escape radius is 
increased up to Rs = 100 (Fig. 4) the fractal dimension value 
is the same as for an escape radius Rs = 10 (Fig. 3). Thus 



Figure 4: Number of boxes Nt covering at scale b non-escaped orbits structure 
(strange repeller) appearing in (x, y) - Poincare section of the Jacobi Hamilto¬ 
nian Box-counting computation is performed as in Fig. 3. Box-counting 
results are shown for Poincare section of orbits staying in the r < Rs = 100 
disk after t = 3 (black crosses), t = 10 (red squares), t = 30 (green circles), 
r = 50 (blue triangles). Keeping orbits staying in the r < Rs = 100 disk 
after time t = 3, 10,30,50 we obtain a strange repeller fractal dimension (a) 
D = 1.8908 + 0.00876,1.8900 + 0.00858,1.8874 + 0.00754,1.8799 + 0.00480 
for// = 0.3 and C = 3 (see Fig. la,b,c), (b) D = 1.8916 + 0.0129,1.8911 + 
0.0127,1.8869 + 0.0110,1.8786 + 0.0086 for// = 0.4 and C = 3 (see Fig. lOb). 
Fits have been performed in the scale range 2^ < l/Z? < 2^. We used A = 10^ 
test particles initially distributed at random in the 1.3 < r < 2.5 ring. The fractal 
dimension has been computed with (a) Nt >3 = 38090345, A^>io = 18470667, 
A, >30 = 7588914, A, >50 = 4574705 points, (b) A, >3 = 27206778, A, >10 = 
14259496, A, >30 = 6542192, A, >50 = 4286763 points. 


the obtained value of D is not affected by the escape cut-off 
distance Rs. Also, as seen in Fig. 4, the fractal dimension re¬ 
mains practically the same if we consider strange repeller struc¬ 
tures obtained after time f = 3, 10, 30 and 50. Hence even for 
short times the strange repeller structure is already well defined 
and perdures for greater times since D is constant (at least here 
up to f = 50). Throughout this work, for the sake of clarity 
we choose to present Poincare sections for orbits staying in the 
r < Rs = 10 disk after time t = 10. 

Fig. 5 shows (x, y) - Poincare sections for the mass param¬ 
eter A = 0.3 and for different Jacobi constants from C = 2.6 
to C = 4. The strange repeller structure constituted by non¬ 
escaping orbits is progressively expelled from the primaries 
vicinity as C increases. At C = 2.6, 3, 3.2, 3.4 (Fig. 5a, b, c,d) 
non-escaping trajectories may still pass close by each one of the 
primaries. The strange repeller still dominates the phase por¬ 
trait with a fractal dimension decreasing down to Z) 1.84 for 
C = 3.4 (Fig. 3J). For greater values C = 3.6, 4 (Fig. 5e, f) the 
forbidden zone insulates the immediate vicinity of the primaries 
from trajectories coming from regions beyond r ~ 1. Regular 
and chaotic trajectories corresponding to particles gravitating 
one or the two primaries are confined in the very central region 
( Valtonen and Karttunen[ [2006] ). The strange repeller is con¬ 
fined in a narrow ring located in the region r ~ 1.5. Beyond 
that region we observe nearly stable circular orbits (blue dots) 
corresponding to particles gravitating the whole binary with a 
radius r ~ 2. For these orbits stable means that these orbits have 
not escaped from the disk r < Rs during the whole integration 
duration ts. 

Unless otherwise stated, we have deliberately omitted the 
class of orbits crashing primaries. According to |Nagler| ( |2004 


2005) the crash basin scales as a power law where is the 
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Figure 6: (x,y) - Poincare sections of the Jacobi Hamiltonian ^ with r = 0 
and 0 < 0 for primary bodies with radius = 0.01 and for // = 0.5 and C = 1. 
The panels show traces of orbits associated with particles escaping at time (a) 
^esc ^ 0 (^i.6. all computed orbits are shown), (Z?) > 0.01, (c) > l,(d) 

tesc >10. The gray scale bar shows the time when particles pass through the 
Poincare section. Light gray (dark gray) points have been obtained at r ^ 1 
(t ^ 10^). Red (blue) points have been obtained from orbits crashing on the ju 
(l-ju) primary mass. Initially 10^ particles have been randomly distributed in 
the 1.3 < r < 2.5 ring. Fig. 5 (4th panel for = 0.01) in |Nagler|j2004) is 
similar to panel (a). 


Figure 5: (x,y) - Poincare sections of the Jacobi Hamiltonian with r = 0 
and 0 < 0 for // = 0.3 and Jacobi constant (a) C = 2.6, (b) C = 3, (c) C = 3.2, 
(d) C = 3.4, (e) C = 3.6, and (/) C = 4. Red regions are forbidden since 
there x^ + y^ < 0. Black dots represent non-escaped orbits staying inside the 
r < = 10 region after time t = 10. Blue dots represent bounded orbits. 

The red (blue) star * (*) gives the position of the primary (secondary) mass as 
in Fig. 1. Poincare sections have been obtained analyzing orbits of N = 
(a,b,c) and N = 10^ {d,e,f) particles initially placed in the 1.3 < r < 2.5 
region. Particles are considered as escaped once r > Rs. 


radius of the primary mass and a ~ 0.5. In this work we choose 
a radius of rt, = 10“^ for the two primaries which gives two 
percent of crashing orbits and an area of about one percent for 
crash basin not visible in the presented Poincare sections. The 
sets of non-escaping orbits shown in Fig. 1 are also distinguish¬ 
able in the studies ( |Nagler| |2004[ |2005| ) devoted to crashing 
orbits but not studied in details. For example Fig. Id presents 
a Poincare section for the Copenhagen problem case (// = 0.5) 
with C = 3 which is similar to the Poincare section presented 
in the study (|Nagler| |2004|) Fig. 3 right column middle row for 
C = 2.85. 


The time evolution of density of non-escaped particles is 
shown in Fig. 6 for yu = 0.5 and Fig. 7 for yu = 0.3 for the case of 
primary bodies of relatively large radius rt, = 0.01 (such a size 
is also available in Fig. 5 in |Nagler| ( |2004| )). These data clearly 


show that the strange repeller structure is established on rather 
short time scales with f ~ 1. We also see that for such a value 
of rij the measure of crashed orbits gives a visible contribution 
to the measure of non-escaping orbits of strange repeller. 

Finally we note that in our computations we determined the 
fractal dimension D of trajectories non-escaping in future (see 
a similar situation considered for the Chirikov standard map 
with absorption in Ermann, L. and Shepelyansky, D. L.] 2Q10| ). 
According to the general relations known for the strange frac¬ 
tal sets in dynamical systems the fractal dimension of the 
invariant repeller set (orbits non-escaping neither in the future 
nor in the past) satisfies the relation Dq = 2(D - 1) ( |Qtt[|1993) . 
Thus for the typical value we have in Fig. 1 at yu = 0.3 with 
Z) 1.87 we obtain Dq ~ 1.74. 

We can expect that the spiral fractal structure, clearly present 
in the plane (x,y) of the Poincare sections (see e.g. Fig. 1), will 
give somewhat similar traces for the surface (or area) density 
Ps = dNtjdxdy of particles Nt remaining in the system at an 
instant moment of time t. A typical example of surface density, 
corresponding to Fig. \a,b,c, is shown in Fig. 8. Indeed, we 
find a clear spiral structure of ps{x, y) similar to the spiral struc¬ 
ture of the strange repeller of Fig. \a, h, c. Of course, here in 
Fig. 8 we have the projected density of particles in (x,y) plane 
taken at all values of r (and not only at r = 0 as in Fig. 1), This 
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Figure 7: (x, y) - Poincare sections of the Jacobi Hamiltonian with r = 0 and 
0 < 0 for primary bodies with radius = 0.01 and for // = 0.3 and C = 3. The 
panels show (a) all points after t = 0, (b) all points after t = 0.01, (c) all points 
after r = I, id) all points after r = 10. The gray scale bar shows the time when 
particles pass through the Poincare section. Light gray (dark gray) points have 
been obtained at r ^ 1 (f ^ 10^). Red (blue) points have been obtained from 
orbits crashing on the ji (l-n) primary mass. Initially 10^ particles have been 
randomly distributed in the 1.3 < r < 2.5 ring. 


Figure 8: {a) Snapshot of the remaining particles at r = 10 for ju = 0.3 and 
C = 3. The number of particles initially injected in the 1.3 < r < 2.5 ring is 
N = 10^, at time r = 10 there are Nt = 13302225 particles remaining inside 
the circle r < Rs = 10, colors show the surface density of particles ps in the 
plane (x,y), color bar gives the logarithmic color scale of density with levels 
corresponding to a proportion of the Nt remaining particles; (b) same as panel 
(a) but on a smaller scale; (c) linear density p(r) = dNtIdr, (d) angle averaged 
surface density p(r)lr = 2n < ps(r) >, the dashed line shows the slope oc 1/r^/^. 


leads to a smoothing of the fractal structure of the Poincare sec¬ 
tion but the spiral distribution of density ps remains well visi¬ 
ble. We also note that the angle averaged density of remaining 
particles < Ps(r) >oc has a radial dependence on r be¬ 

ing similar to those found for the dark matter density obtained 
in the symplectic simulations of scattering and capture of dark 


matter particles on binary system (see e.g. Fig. 4a in Lages and 
Shepelyansky[|2Q13}|Rollin et aH 2015b ). 


4. Poincare recurrences and probability decay 

We determine numerically the probability P(t) to stay inside 
the system for time larger than t. For that N particles are ran¬ 
domly distributed in the range 1.3 < r < 2.5 and then the sur¬ 
vival probability P(t) is defined as the ratio P(t) = Nt/N, where 
Nt is the number of particles remained inside the system with 
r < Rs = 10 at times larger than t. This survival probability 
is proportional to the integrated probability of Poincare recur¬ 
rences (see e.g. [Chirikov and Shepelyanskyj |1999t |Frahm, K. 
M. and Shepelyansky, D. L.[]2010 ). 

The typical examples of the decay P(t) are shown in Fig. 9. 
At relatively small time scales with f < 100 the decay can be 
approximately described by an exponential decay with a de¬ 
cay time Ts ~ 10. Indeed, for the strange dynamical sets (e.g. 
strange attractors) one obtains usually an exponential decay 


since there are no specific sticking regions in such strange sets 
( Chirikov and Shepelyansky] |1984| ). 

However, at larger time scales t > 100 we see the appearance 
of the algebraic decay of probability corresponding to the alge¬ 
braic statistics of Poincare recurrences discussed for symplec¬ 
tic maps (see e.g. Chirikov and Shepelyansky 1999[ Cnstadorol 


and Ketzmerick[ 20081 Frahm, K. M. and Shepelyansky, D. L. 


2010). When the decay time is measured in number of cross 
ings of the Poincare section n we obtain the Poincare exponent 
P of this decay p = 1.49 being close to the values p ^ 5 found 

in the symplectic maps. However, if the time is measured in 
number of rotations of binaries t/lTi then we find a somewhat 
large value of p (see Fig. 9). We explain this deviation a bit 
later. 

The properties of orbits surviving in the system for long 
times are shown in Fig. 10. We see that such orbits are concen¬ 
trated in the vicinity of critical KAM curves which separate the 
orbits of strange repeller from the integrable islands with KAM 
curves. This is exactly the situation discussed in the symplec¬ 
tic maps (see e.g. Chirikov and Shepelyansky 1999[ Cnstadorol 


and Ketzmerick 20081 Frahm, K. M. and Shepelyansky, D. L. 


2010| ). Thus the asymptotic decay of survival probability is de¬ 
termined by long time sticking orbits in the vicinity of critical 
KAM curves. The detail analytical explanation of this generic 
phenomenon is still under debates (see e.g. |M eiss l |20T5l l. 
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Figure 9: Survival probability P of particles as a function of time t (left panels, 
binary period is 2n) and as a function of the number n of successive Poincare 
section crossings (right panels) for C = 3 and binaries with mass parameter 
jj. = 0.3 (black curve), // = 0.4 (red curve), id = 0.5 (blue curve). Survival 
probabilities are shown in semi-log scale (top panels) and in log-log scale (bot¬ 
tom panels). Dashed lines show {a) exponential decay P oc exp(-r/T5) with 
l/r^ = 0.13, (b) algebraic decay Pit) oc 1//^ with the Poincare exponent 
yS = 1.82, (c) exponential decay Pin) oc exp(-n/T5) with l/r^ = 0.07, id) 
algebraic decay Pin) oc 1/^^ with/5 = 1.49. Initially 10^ particles have been 
randomly distributed in the 1.3 < r < 2.5 ring. Escape radius is Rs = 10. 


Figure 11: Survival probability P of particles as a function of time t (left 
panels, binary period is 2n) and as a function of the number n of successive 
Poincare section crossings (right panels) for n = 0.3 and C = 2.6 (black curve), 
C = 3 (red curve), C = 3.4 (blue curve). Survival probabilities are shown in 
lin-log scale (top panels) and in log-log scale (bottom panels). Blue dashed 
lines show exponential decay Pit) oc exp(-r/T;y) with l/r^ = 0.04 (a,b) and 
Pin) oc Qxpi-nlTg) with l/r^ = 0.03 (c,d). Gray dashed lines show ia) expo¬ 
nential decay Pit) with l/r^ = 0.14, ib) algebraic decay Pit) with /5 = 1.87, (c) 
exponential decay Pin) with 1 /r^ = 0.09, id) algebraic decay Pin) with/5 = 1.6. 
Initially 10^ particles have been randomly distributed in the 1.3 < r < 2.5 ring. 
Escape radius is Rs = 10. 


At fixed ju = 03 the variation of decay properties of P(t) 
with the Jacobi constant C is shown in Fig. 11. For C = 3.4 
the exponential decay of P(t) is much slower than for the cases 
C = 3.6 and C = 3 (there is a factor 3.5 between corresponding 
characteristic time scales r^). We attribute this behavior to the 



Figure 10: (x,y) - Poincare sections of the Jacobi Hamiltonian with r = 0 
and 0 < 0 for C = 3 in the case of ia) an = 0.3 binary and ib) an = 0.4 binary. 
Red regions are forbidden since there < 0. Black dots represent non- 

escaped orbits staying inside the r < = 10 region after time r = 10. Blue 

dots represent bounded orbits on KAM curves inside integrable islands. Green 
plain circles mark out non-escaped orbits remaining inside the disk r < Rg 
after time t = 500. The red (blue) star * (*) gives the position of the primary 
(secondary) mass. Each of Poincare sections is obtained from orbits of N = 
10^ particles initially placed in the 1.3 < r < 2.5 region; these particles are 
considered as escaped once r > Rg. 


fact that the forbidden zone encloses almost completely binary 
components letting for particles only a small route to binary 
components in comparison with cases C = 2.6 and C = 3 for 
which the route to gravitational perturbers is not constrained. 
Meanwhile, particles have a limited access to the chaotic com¬ 
ponent (see Fig. 5d where it forms a narrow peanut shell around 
central stability islands) so that the sticking in the vicinity of 
critical KAM curves is reduced emphasizing the exponential 
decay. For C = 2.6,3 at long time scales (Fig. lid) we have 
a well visible algebraic decay P(n) with yS ~ 1.5 (we may as¬ 
sume that increasing sufficiently the number of initially injected 
particle for the case C = 3.4, blue line in Fig. lid would also 
exhibit at large time scales the same yS ~ 1.5 algebraic decay). 
However, algebraic decay of P(t) for C = 2.6,3 at large time 
scale (Fig. 11b) still have somewhat different value of (3. 

The origin of this difference for P{t) becomes clear from 
Fig. 12 where we show the data similar to those of Figs. 9,11 
at yu = 0.3 but with the escape radius Rg = 100. We see that 
the decay properties of P{n) remain practically unchanged that 
confirms the generic features of obtained results for fd (indeed, 
stability islands do not affect dynamics at r ~ = 100). How¬ 

ever, the value of jd for Pit) is significantly reduced fofd ^ 1.1. 
We explain this by the fact that in the usual time units the mea¬ 
sure of chaotic component at large distances becomes dominant 
and the escape time is determined simply by a Kepler rotation 
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Figure 12: Survival probability P of particles as a function of time t (left panels, 
binary period is 27r) and as a function of the number n of successive Poincare 
section crossings (right panels) for // = 0.3, C = 3 in the case of an escape 
radius Rs = 100. Dashed curves show (a) exponential decay Pit) with l/r^ = 
0.012, ib) algebraic decay Pit) with p = 1.11, (c) exponential decay Pin) with 
l/r^ = 0.062, id) algebraic decay Pin) with = 1.40. Initially 10^ particles 
have been randomly distributed in the 1.3 < r < 2.5 ring. 



10'^ 




10'^ 


period which becomes larger for large r values. This leads to 
the decay exponent yS = 2/3 for Pit) as discussed in |Borgonovi 


et al.| ( |1988| ) and explains the variation of p with Rs. How¬ 


ever, when the time is counted in the number of orbital periods 
the divergence of the orbital period at large r values (or small 
coupling energies) does not affect the decay and we obtain the 
Poincare exponent p ^ 1.5 being independent of Rs. 

We note that the recent studies of survival probability decay 
in the PCRTBP also report the value of p ^ 1.5 ( [Kovacs and 
Regaly 2015| ). 


5. Symplectic map description 


Figure 13: Density of particles in the steady state, (a) Number N of particles in 
the r < 100 region of the binary as a function of time t. We have continuously 
injected N = 10^ particles placed randomly on r = 100 circle with close to 
parabolic but hyperbolic trajectories in the 2-body paradigm. Steady state is 
achieved around t ~ tQ = 10^. The surface density of particles Psix,y) is 
shown for Jacobi constants C = 2.6 (b), C = 3A (c) and C = 4 (e). Each 
surface density has been obtained from 1000 surface density snapshots taken at 
regular time interval after fo = 10^ once steady state is attained. Color bar gives 
the logarithmic color scale of density with levels corresponding to the ratio of 
Nc/Nio where Nc is the number of particles in a given cell of size 0.1 x 0.1 
and Nio is the number of particles in the r < 10 region. Panel (d) shows the 
normalized steady state distribution of Jacobi constants of particles present in 
the r < 100 region. Panel (f) shows the steady state surface density build with 
all particles in r < 100 region. 


Finally we discuss the case when particles in the sidereal ref¬ 
erence frame scatter on the binary with relatively large values 
of perihelion distance ^ » 1. Such a case corresponds to large 
|C| » 1. For ^ » 1 the dynamics of particle in the field of 
two binaries is approximately described by the symplectic map 
of the form w = w Fix); x = x where w = -2E is 


the particle energy, v is the phase of binary rotation (in units of 
In) when the particle is located at its perihelion and Fix) oc fi 
is a periodic fu nction of v ([Petrosky 1986 Chirikov and Vech- 


|eslavov||1989t|^evchenko||201H|Rollin et al.||2015a|b| ); bars 

above w and v mark new values after one rotation around the 
binary. The amplitude F^ax decreases exponentially with in¬ 
creasing q. Usually, one considers the case of yu 1 (e.g. Sun 
and Jupiter) but our studies show that this form remains valid 
even for ~ 0.5 if ^ » 1. At /i 1 we have F^ax ^ 1 


and the escape time ti becomes very large ti oc 1 IF^ax^ being 
much larger than the Lyapunov time scale. In this situation the 
fractal dimensions D and Dq are very close to D = Dq = 2 
( Qtt||1993]|Ermann, L. and Shepelyansky, D. L.[|201Q| ) and the 
fractal effects practically disappear. Due to that this case is not 
interesting for the fractal analysis. 


6. Discussion 

We analyzed the PCRTBP dynamical system and showed that 
for moderate mass ratio of primary bodies ju - 0.5 the Poincare 
section is characterized by a strange repeller with the fractal 
dimension having typical values D ^ 1.87 (Dq ~ 1.74) at mod¬ 
erate values of the Jacobi constant C ~ 2. At the same time cer¬ 
tain islands of integrable motion are still present. Such a struc- 
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ture of the Poincare section leads to an exponential decay of 
survival probability in the system on short time scales followed 
by the algebraic decay with the Poincare exponent (3 ^ 5 be¬ 

ing similar to the values known for the statistics of Poincare 
recurrences in generic symplectic maps. For the small mass ra¬ 
tio yu ^ 1 the escape times becomes very large and the fractal 
dimension becomes close to the usual value D = 2. 


It is interesting to note that the strange repeller structure (see 
Figs. 1,5,6,7,10) reminds the structure of spiral galaxies (see 


e.g. Milky Way representation at European Southern Observa- 
tory websit^2013] r. Intact the fractal features of galaxies have 


been studied extensively by various groups ( [Elmegreen and 
Elmegreen] [2001 [ [Sanchez et al.| [2007 1 12010| ) which obtained 


from observation data the fractal dimensions for the plane den¬ 
sity being around Dg ^ 1.3 - 1.7 for different galaxies. Some of 
these values (e.g. [Sanchez et al.||2007| with Dg ^ 1.7) are sim¬ 


ilar to those obtained here for the Poincare section of PCRTBP. 
Our results show that the spiral fractal structure of the Poincare 
section of PCRTBP leads to a spiral structure of global den¬ 
sity of particles ps remaining in the system (see Pig. 8). Thus 
we make a conjecture that the spiral structure of certain galax¬ 
ies can be linked to the underlying spiral fractals appearing 
in the dynamics of particles in binary systems. Of course, 
such a conjecture requires more detailed analysis beyond the 
present paper such as taking into account the third dimension. 
Also whether or not the spiral structure remains in the steady 
state is an important question. To bring some elements of re¬ 
sponse to the latter question we have injected at r = 100 par¬ 
ticles with close-to-parabolic hyperbolic trajectories in 2-body 
paradigm. We use a Maxwellian distribution for initial veloci¬ 
ties f(v)dv = ^J54/7^ exp dv with u = 0.035 

and a homogeneous distribution for perihelion q g [0; 3] pa¬ 
rameters. We checked that our results in steady state are in¬ 
dependent of initial distributions. We present on Pig. 13 the 
surface density of particles p^(v,y) in the steady state attained 
around ~ 10^ (Pig. 13a). We clearly see that the two arms 
spiral structure remains in the steady state, we show examples 
for particles with trajectories characterized by (b) C = 2.6, (c) 
C = 3.4, (e) C = 4 (example not shown for C = 3 is analogous 
to Pig 8a). However when all injected particles are considered 
independently of the Jacobi constant C we obtain a blurred spi¬ 
ral (Pig. 13f). The contributions of different Jacobi constants 
are consequently mixed according to the distribution Pig. 13d. 
We remark that the results presented in Pig. 13 share similarities 
with those obtained in the study of chaotic trajectories in spi¬ 
ral galaxies ( [Harsoula et al.[ [20111 [Contopoulos and Harsoula 


[2012[[MT^ . On the basis of obtained results we make a con¬ 
jecture about existence of certain links between observed frac¬ 
tal dimensions of galaxies ( [Elmegreen and Elmegreenj [2001 [ 
Sanchez et al. 2007 [ [2010[ ) and fractal spiral repeller structure 


studied here. We think that the further studies of fractal struc¬ 
tures in binary systems will bring new interesting results. 
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